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We investigate tiie iiydrostatic equilibrium of stellar structure by taking into account the modi- 
fied Lane-Emden equation coming out from /(_R)-gravity. Such an equation is obtained in metric 
approach by considering the Newtonian limit of /(i?)-gravity, which gives rise to a modified Poisson 
equation, and then introducing a relation between pressure and density with polytropic index n. The 
modified equation results an integro-differential equation, which, in the hmit f{R) — >■ R, becomes 
' the standard Lane-Emden equation. We find the radial profiles of gravitational potential by solving 

^SJ ' for some values of n. The comparison of solutions with those coming from General Relativity shows 

^ , that they are compatible and physically relevant. 
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I. INTRODUCTION 



Extended Theories of Gravity (ETG) [1| are a new paradigm of modern physics aimed to address several shortcom- 
ings coming out in the study of gravitational interaction at ultra-violet and infra-red scales. In particular, instead 
of introducing unknovifn fluids, the approach consists in extending General Relativity (GR) by taking into account 
generic functions of curvature invariants. These functions can be physically motivated and capable of addressing 
phenomenology at galactic, extragalactic, and cosmological scales [5]. 
^ . This viewpoint does not require to find out candidates for dark energy and dark matter at fundamental level (not 
' detected up to now), but takes into account only the observed ingredients (i.e. gravity, radiation and baryonic matter), 
changing the l.h.s. of the field equations. Despite of this modification, it is in agreement with the spirit of GR since 
the only request is that the Hilbert-Einstein action should be generalized asking for a gravitational interaction acting, 
[ in principle, in different ways at different scales but preserving the robust results of GR at local and Solar System 
scales (see [l| for a detailed discussion). This is the case of /(i?)-gravity which reduces to GR as soon as f{R) — >■ R. 

Other issues as, for example, the observed Pioneer anomaly problem can be framed into the same approach Q 
and then, apart the cosmological dynamics, a systematic analysis of such theories urges at short scales and in the low 
energy limit. 

On the other hand, the strong gravity regime [sj is another way to check the viability of these theories. In 
general the formation and the evolution of stars can be considered suitable test-beds for Alternative Theories of 
^ . . Gravity. Considering the case of /(i?)-gravity, divergences stemming from the functional form of /(i?) may prevent 
' the existence of relativistic stars in these theories , but thanks to the chameleon mechanism, introduced by Khoury 
and Weltman [?[, the possible problems jeopardizing the existence of these objects may be avoided 3. Furthermore, 
there are also numerical solutions corresponding to static star configurations with strong gravitational fields 9] where 
the choice of the equation of state is crucial for the existence of solutions. 

Furthermore some observed stellar systems are incompatible with the standard models of stellar structure. We refer 
to anomalous neutron stars, the so called "magnetars" [l^ with masses larger than their expected Volkoff mass. It 
seems that, on particular length scales, the gravitational force is larger or smaller than the corresponding GR value. 
For example, a modification of the Hilbert-Einstein Lagrangian, consisting of terms, enables a major attraction 
while a RafsR"^ term gives a repulsive contribution [ll|. Understanding on which scales the modifications to GR are 
working or what is the weight of corrections to gravitational potential is a crucial point that could confirm or rule 
out these extended approaches to gravitational interaction. 
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The plan of paper is the fohowing: In Sec. HIl we review briefly the classical hydrostatic problem for stellar 
structures. In Sec. IIIII we derive the Newtonian limit of /(i?)-gravity obtaining the modified Poisson equation. The 
modified Lane-Emden equation is obtained in Sec. HV] and its structure is compared with respect to the standard 
one. In Sec|Vl we show the analytical solutions of standard Lane-Emden equation and compare them with those 
obtained perturbatively from /(i?)-gravity. With help of plot we can compare between them all results. Discussion 
and conclusions are drawn in Sec. I VII . 

II. HYDROSTATIC EQUILIBRIUM OF STELLAR STRUCTURES 

The condition of hydrostatic equilibrium for stellar structures in Newtonian dynamics is achieved by considering 
the equation 



dp d$ 

where p is the pressure, — $ is the gravitational potential, and p is the density [l^ . Together with the above equation, 
the Poisson equation 



gives the gravitational potential as solution for a given matter density p. Since we are taking into account only static 
and stationary situations, here we consider only time-independent solutions ^ . In general, the temperature r appears 
in Eqs. ([T]) and ^ the density satisfies an equation of state of the form p = p{p,T). In any case, we assume that 
there exists a polytropic relation between p and p of the form 



p — Kp"* 



(3) 



where K and 7 are constant. Note that $ > in the interior of the model since we define the gravitational potential 
as — The polytropic constant K is fixed and can be obtained as a combination of fundamental constants. However 
there are several realistic cases where K is not fixed and another equation for its evolution is needed. The constant 
7 is the polytropic exponent . Inserting the polytropic equation of state into Eq. ([T|) , we obtain 



— = 7XpT — 
dr dr 



For 7 7^ 1, the above equation can be integrated giving 



(4) 



7^ 



.7-1 



P 



(5) 



where we have chosen the integration constant to give $ = at surface (p = 0). The constant n is called the polytropic 
index and is defined as 71 = . Inserting the relation ([S]) into the Poisson equation, we obtain a differential equation 
for the gravitational potential 



dr^ 



r dr 



(6) 



Let us define now the dimensionless variables 



(7) 



^ The radius r is assumed as the spatial coordinate. It varies from r 



= at the center to r = ^ at the surface of the star 
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where the subscript c refers to the center of the star and the relation between p and $ is given by Eq. ([5|). At the 
center (r = 0), we have 2 = 0,$ = $cj P = Pc and therefore w — 1. Then Eq. © can be written 



(fw 2 dw 

dz^ z dz 



This is the standard Lane-Embden equation describing the hydrostatic equihbrium of stellar structures in the Newto- 



nian theory 12 | 



III. THE NEWTONIAN LIMIT OF /(7?)-GRAVITY 

Let us start with a general class of Extended Theories of Gravity given by the action 



A = j d^x^[f(R) + XC„,] (9) 

where f{R) is an analytic function of the curvature invariant R. Cm is the minimally coupled ordinary matter 
Lagrangian density. In the metric approach, the field equations are obtained by varying the action (|9]) with respect 
to g^v. We get 

/ ^g^u f-.^u ^" Q^y^I — ^ T^jy 

(10) 

3n/' + /'i?- 2/ ^ XT 



where the second equation is the trace of the field equations. Here, T^ji/ = ^ ^ ^^^^ ™ is the the energy- 



momentum tensor of matter; T ~ T" ^ is the trace; /' — □ = lo-''^ the d'Alembert operator and X — 8t:G. 

We assume c = 1 is adopted. The conventions for Ricci's tensor are = R" ^cru] the Riemann tensor is R"' p^u = 
^ + .... The affine connections are the ChristofFel's symbols of the metric F^^^ — \g^'^ {gacr,ii + gi3a,a — ffa^.o-)- The 
signature is (H ). 

In order to achieve the Newtonian limit of the theory the metric tensor have to be approximated as follows 

/l-2$(t,x)+0(4) 0(3) \ 

9^.^- \ (11) 
\ 0(3) -%+0(2)y 

where 0{n) (with n — integer) denotes the order of the expansion (see [l^ for details). The set of coordinates^ 
adopted is = {t,x^ ,x'^ ,x^). The Ricci scalar formally becomes 

R - i?(2)(i,x) +0(4) (12) 
The n-th derivative of Ricci function can be developed as 



r{R) - r(i?(2) + 0(4)) ^ /«(o) + /"+i(0)i?(2) + 0(4) (13) 

here i?*^"^ denotes a quantity of order 0{n). From lowest order of field Eqs. ([TU)) . we have /(O) = which trivially 
follows from the above assumption ([TT]) for the metric. This means that the space-time is asymptotically Minkowskian 
and we are discarding a cosmological constant term in this analysis'^. Eqs. (|10p at C'(2)-order, that is at Newtonian 
level, are 



^ The Greek index runs between and 3; the Latin index between 1 and 3. 

^ This assumption is quite natural since the contribution of a cosmological constant term is irrelevant at stellar level. 
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i?(2) -^B^- /"(0)Ai?(2) = XT^l 



(0) 

(14) 

3/"(0)Ai?(2) _ 7^(2) ^ ;fr(0) 



where A is the Laplacian in the flat space, — — A$(t,x) and, for the sake of simplicity, we set /'(O) — 1. We 
recall that the energy-momentum tensor for a perfect fluid is 

T^_,u = {e + p)Uf_,u^-pg^^ (15) 

where p is the pressure and e is the energy density. Being the pressure contribution negligible in the field equations 
in Newtonian approximation, we have 

A$ + ^ + /"(0)Ai?(2) = -A-p 

(16) 

3/"(0)Ai?(2) + i^(2) ^ _p^p 

where /9 is now the mass density^ . We note that for /"(O) = we have the standard Poisson equation: A<i> = —AnGp. 
This means that as soon as the second derivative of /(i?) is different from zero, deviations from the Newtonian limit 
of GR emerge. 

The gravitational potential — $, solution of Eqs. ([T6|. has in general a Yukawa-like behavior depending on a 
characteristic length on which it evolves .133- Then as it is evident the Gauss theorem is not valid'"^ since the force 
law is not cx |x|~2. The equivalence between a spherically symmetric distribution and point-like distribution is not 
valid and how the matter is distributed in the space is very important [isl - fisj . 

Besides the Birkhoff theorem results modified at Newtonian level: the solution can be only factorized by a space- 
depending function and an arbitrary time-depending function [Tsj . Furthermore the correction to the gravitational 
potential is depending on the only first two derivatives of f{R) in R — 0. This means that different analytical 
theories, from the third derivative perturbation terms on, admit the same Newtonian limit (isl [l3|. 

Eqs. ([TCI) can be considered the modified Poisson equation for /(i?)-gravity. They do not depend on gauge condition 
choice |15| . 

IV. STELLAR HYDROSTATIC EQUILIBRIUM IN /(i?)-GRAVITY 

From the Bianchi identity, satisfied by the field Eqs. ((TU)) ). we have 

If the dependence on the temperature r is negligible, i.e. p — p{p), this relation can be introduced into Eqs. (jl6p . 
which become a system of three equations for p, $ and i?*-^'' and can be solved without the other structure equations. 

Let us suppose that matter satisfies still a polytropic equation p — K p^ . If we introduce Eq.(l5]) into Eqs. ([T6l) we 
obtain an integro-differential equation for the gravitational potential — <I>, that is 

A$(x) ^^$(x)" = / dVg(x,x')$(x')" (18) 

3 6 J 

2 g-m|x-x' 

where tj(x, x') = is the Green function of the field operator Ax — for systems with spherical 

A-K |x — x'l 

symmetry and — (for details see [13, [l5[). The integro-differential nature of Eq. p^ is the proof of the 

3/ (0) 

non- viability of Gauss theorem for /(i?)-gravity. Adopting again the dimensionless variables 



* Generally it is e = pc^. 

^ It is worth noticing that also if the Gauss theorem does not hold, the Bianchi identities are always valid so the conservation laws are 
guaranteed. 
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where 



M 

(19) 



is a characteristic length hnked to stellar radius ^, Eq. (fT8| becomes 

^+2d^^^^^^„ ^ /'^'"dz'z'(e-«ol— 'l-e-"«ol^+^'llz.(zr (21) 

(ix^ z dz 8 z I J 

which is the modified Lane-Emden equation deduced from /(i?)-gravity. Clearly the particular /(i?)-model is specified 
by the parameters m and ^o- If m — >■ oo {i.e. f{R) — > R), Eq. (I^TI) becomes Eq. (|H]). We are only interested in 
solutions of Eq. (1211) that are finite at the center, that is for z = 0. Since the center must be an equilibrium point, 
the gravitational acceleration |g| = —d^/dr oc dw/dz must vanish for ^'(0) — 0. Let us assume we have solutions 
w{z) of Ea.(pT|) that fulfill the boundary conditions w(0) = 1 and w(^/^o) = 0; then according to the choice (|19p . 
the radial distribution of density is given by 

p(|x|) = p^w" , pe = A„$,« (22) 

and the pressure by 

p(|x|) = , = i^pj (23) 

For 7 = 1 (or n — oo) the integro-differential Eq. (1211) is not correct. This means that the theory does not contain 
the case of isothermal sphere of ideal gas. In this case, the polytropic relation is p = K p. Putting this relation into 
Eq.dni) we have 

^ = Inp-lnp, p = Pce"^"" (24) 

where the constant of integration is chosen in such a way that the gravitational potential is zero at the center. If we 
introduce Eq. (p4)) into Eqs. (|16p . we have 

A$(x) + ^ rr^Xp^ f ^3^g^^^^>yH^')/K (25) 



Assuming the dimensionless variables z = and w{z) = where = ^ -iXp ' ^"^^ f25)) becomes 

d^^;(z) ^ 2 rfu;(z) ^ ^ m^il f^/^^ ^, f _™^,,,_,'| _ 1 ^(,') ^^e) 

dz^ z dz 8 z Jq \ J 

which is the modified "isothermal" Lane-Emden equation derived /(i?)-gravity. 



V. SOLUTIONS OF THE STANDARD AND MODIFIED LANE-EMDEN EQUATIONS 

The task is now to solve the modified Lane-Emden equation and compare its solutions to those of standard New- 
tonian theory. Only for three values of n, the solutions of Eq.(I5]) have analytical expressions [l^] 
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^ 4")^(,) = i-_ 



/ N sinz . _. 

1 ^ «;g],(z) = — (27) 

1 



n = 5 — > Wq]j(z) = 



1 + — 



We label these solution with on since they agree with the Newtonian limit of GR. The surface of the polytrope of 
index n is defined by the value z = z^^\ where p = and thus w = 0. For n = and n = 1 the surface is reached 
for a finite value of z'^"'\ The case n = 5 yields a model of infinite radius. It can be shown that for rt < 5 the radius 
of polytropic models is finite; for n > 5 they have infinite radius. From Eqs. (|27p one finds z^jj, = a/G, '^g\i ^ ^ ^'^'^ 
ZqI^ = oo. A general property of the solutions is that z^"^ grows monotonically with the polytropic index n. In Fig. 

[T]we show the behavior of solutions w|^^ for n = 0, 1, 5. Apart from the three cases where analytic solutions are 
known, the classical Lane-Emden Eq. ^ has to be be solved numerically, considering with the expression 



41w = E«^^-' (28) 



for the neighborhood of the center. Inserting Eq. ipS)) into Eq. and by comparing coefficients one finds, at lowest 
orders, a classification of solutions by the index n, that is 

+ + (29) 

The case 7 = 5/3 and n = 3/2 is the non-relativistic limit while the case 7 = 4/3 and n = 3 is the relativistic limit 
of a completely degenerate gas. 

Also for modified Lane-Emden Eq. (|2ip . we have an exact solution for n = 0. In fact, it is straightforward to find 
out 



sinh m^oz 



m£,oz 



(30) 



where the boundary conditions w{0) = 1 and w'{0) = are satisfied. A comment on the GR limit (that is 
f{R) — > R) of solution (1301) is necessary. In fact when we perform the limit m oo, we do not recover exactly 

Wq]^{z). The difference is in the definition of quantity ^o- In /(i?)-gravity we have the definition ((20)) while in GR it 
I 2 

is ^0 = \ ^TT' since in the first equation of (|16p. when we perform f{R) — > R, we have to eliminate the trace 

y XAn^c 

equation condition. In general, this means that the Newtonian limit and the Eddington parameterization of different 
relativistic theories of gravity cannot coincide with those of GR (see [11] for further details on this point). 
The point is calculated by imposing w^^^^^ ^^f(R) ) ~ ^ considering the Taylor expansion 

sinhmfnz l.^s9,„.^s/i , , 

mt^QZ 

we obtain z'-'^-' = — — Since the stellar radius £ is given by definition £ = £q z^'^^ , we obtain the 

V3 + (1 + me)e-'< ^ & y 50 

constraint 



^= / 1 (32) 

2ttG ^ /i ^ i+.Ti; ^_„,f 
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By solving numerically the constraint^ Eq. (l32p . we find the modified expression of the radius ^. If m — c« we have 
the standard expression ~ \[^^ valid for the Newtonian limit of GR. Besides, it is worth noticing that in the 

/(i?)-gravity case, for n = 0, the radius is smaller than in GR. On the other hand, the gravitational potential — $ 
gives rise to a deeper potential well than the corresponding Newtonian potential derived from GR [lij . 
In the case rt = 1, Eq. (PT|) can be recast as follows 



+ w{z) = 



dz' \ e-™«ol^-^'l - ^-mio\^+z'\ 1 ^i^^r^ 



(33) 



where w = zw. If we consider the solution of ([55)) as a small perturbation to the one of GR, we have 

^(1) (z) - w^r^l(z) + e-"'^Aw'^^^ (z) 

f(R)\ J GR\ I f(R)\ > 



(34) 



The coefficient e ™^ < 1 is the parameter with respect to which we perturb Eq. (|33|) . Besides these position ensure 
us that when m — ^ cx) the solution converge to something like w'q^^z'). Substituting Ea. ([M| in Eg. ([551), we have 



d^AzD^i) (z) 
and the solution is easily found 



dz' 



-m4o|z- 



(35) 



smz 



1 + 



(l + m^Co^) 



1 



2e 



1 + m^Co' 



(cosf/^o + ?Ti&sin^/^o) 



8(l + m2eo') 



1 + m^^o 



/ ^ 1^ , . , , sinhm£o2 
— ^ (cos C/t + "^Co sm C/4o ) : V cos z 



m^QZ 



(36) 



Also in this case, for m — ^ cx), we do not recover exactly w[i]j(z). The reason is the same of previous n = case 
[l6j . Analytical solutions for other values of n are not available. 

To conclude this section, we report the gravitational potential profile generated by a spherically symmetric source 
of uniform mass with radius ^. We can impose a mass density of the form 



(37) 



where Q is the Heaviside function and M is the mass [Tj,[Tl]- By solving field Eqs. (ITBl) inside the star and considering 
the boundary conditions w(0) = 1 and w'(0) = 0, we get 



(^) 



1 



2^ m^£, 



2t3 



1 



^0 z^ e ™^(1 + rn^) sinhm^o^ 



2^ m^^-^ 



2^3 



2t3 



miaz 



(38) 



In the limit m — )■ oo, we recover the GR case wgr{z) = 1 — In Fig. [T]we show the behaviors of w^"^^ (z) 



and u'^.J^jj [z) with respect to the corresponding GR cases. Furthermore, we plot the potential generated by a uniform 
spherically symmetric mass distribution in GR and /(i?)-gravity and the case Wg]j(z). 



VI. DISCUSSION AND CONCLUSIONS 



In this paper the hydrostatic equilibrium of a stellar structure in the framework of /(i?)-gravity has been considered. 
The study has been performed starting from the Newtonian limit of /(i?)-field equations. Since the field equations 



In principle, there is a solution for any value of m. 



8 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

1 2 3 4 5 

z 

FIG. 1: Plot of solutions (blue lines) of standard Lane-Emden Eq. ([S]): ^^^^(z) (dotted line) and WqI^{z) (dashed line). The 
green line corresponds to WqI^{z). The red lines are the solutions of modified Lane-Emden Eq. (|2ip : w^^J^^ (z) (dotted line) 
and (z) (dashed line). The blue dashed-dotted line is the potential derived from GR {woRiz)) and the red dashed-dotted 

line the potential derived from /(i?)-gravity (^^^^^(z)) for a uniform spherically symmetric mass distribution. The assumed 
values are = 1 and m^o = .4. From a rapid inspection of these plots, the differences between GR and f{R) gravitational 
potentials are clear and the tendency is that at larger radius z they become more evident. 

satisfy in any case the Bianchi identity, we can use the conservation law of energy-momentum tensor. In particular 
adopting a polytropic equation of state relating the mass density to the pressure, we derive the modified Lane- 
Emden equation and its solutions for ti = 0, 1 which can be compared to the analogous solutions coming from the 
Newtonian limit of GR. When we consider the limit f{R) — ^ R, we obtain the standard hydrostatic equilibrium 
theory coming from GR. A peculiarity of /(i?)-gravity is the non- viability of Gauss theorem and then the modified 
Lane-Emden equation is an integro-differential equation where the mass distribution plays a crucial role. Furthermore 
the correlation between two points in the star is given by a Yukawa-like term of the corresponding Green function. 

These solutions have been matched with those coming from GR and the corresponding density radial profiles have 
been derived. In the case n = 0, we find an exact solution, while, for n — I, we used a perturbative analysis with 
respect to the solution coming from GR. It is possible to demonstrate that density radial profiles coming from f{R)- 
gravity analytic models and close to those coming from GR are compatible. This result rules out some wrong claims 
in the literature stating that /(i?)-gravity is not compatible with self-gravitating systems. Obviously the choice of the 
free parameter of the theory has to be consistent with boundary conditions and then the solutions are parameterized 
by a suitable "wave-length" m = that should be experimentally fixed. 

The next step is to derive self-consistent numerical solutions of modified Lane-Emden equation and build up realistic 
star models where further values of the polytropic index n and other physical parameters, e.g. temperature, opacity, 
transport of energy, are considered. Interesting cases are the non-relativistic limit {n = 3/2) and relativistic limit 
(n = 3) of completely degenerate gas. These models are a challenging task since, up to now, there is no self-consistent, 
final explanation for compact objects (e.g. neutron stars) with masses larger than Volkoff mass, while observational 
evidences widely indicate these objects |1(|. In fact it is plausible that the gravity manifests itself on different 
characteristic lengths and also other contributions in the gravitational potential should be considered for these exotic 
objects. As we have seen above, the gravitational potential well results modified by higher-order corrections in the 
curvature. In particular, it is possible to show that if we put in the action (|9]) other curvature invariants also repulsive 
contributions can emerge [lll.ll5|. These situations have to be seriously taken into account in order to address several 
issues of relativistic astrophysics that seem to be out of the explanation range of the standard theory. 
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